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Abstract 

We study the potential observation at the LHC of CP-violating effects in stop 
production and subsequent cascade decays, gg iiU, ii —>■ txp Xj — > Xi^~^^~^ 
within the Minimal Supersymmetric Standard Model. We study T-odd asymmetries 
based on triple products between the different decay products. There may be a large 
CP asymmetry at the parton level, but there is a significant dilution at the hadronic 
level after integrating over the parton distribution functions. Consequently, even 
for scenarios where large CP intrinsic asymmetries are expected, the measurable 
asymmetry is rather small. High luminosity and precise measurements of masses, 
branching ratios and CP asymmetries may enable measurements of the CP-violating 
parameters in cascade decays at the LHC. 



1 Introduction 



The Minimal Supcrsymmetric Standard Model (MSSM) is a particularly compelling ex- 
tension of the Standard Model, that may soon be explored at the Large Hadron Collider 
(LHC). Current data suggest that, if the MSSM is reahsed in Nature, the supersymme- 
try scale should easily be within reach of the LHC design centre-of-mass energy of 14 
TeV [1,2]. If supersymmetry is discovered, many studies will be required to determine 
the exact details of its realisation. 

The MSSM contains a large number of (as yet) undetermined parameters that may 
have non-zero phases [3,4]. Many of these phases are unphysical in the sense that they can 
be absorbed into the definitions of the fields; however, not all phases can be consistently 
removed in this way. In the neutralino/chargino sector of the complex MSSM, the phase of 
the SU(2) gaugino mass M2 is usually absorbed, whereas the phases of the U(l) gaugino 
mass Ml and the Higgsino mixing parameter n are generally left manifest; this is the 
parameterisation we use. The trilinear couplings Af can also be complex. Studies of 
these CP- violating parameters in sparticle decays and via other properties measurable at 
the LHC will be challenging [5,6]. However, they are extremely important, and provide 
a valuable training ground for exploring the limits of the LHC's capabilities. 

Certain combinations of the CP- violating MSSM phases are constrained by the exper- 
imental upper bounds on the electric dipole moments (EDMs) of the electron, neutron 
and atoms, notably '^^^Tl and -"^^^Hg. Ignoring possible cancellations, the most severely 
constrained individual phase in the MSSM is that of /i, which contributes at the one-loop 
level. For O(IOO) GeV masses, one must require \(f)^\ ^0.1. However, this restriction can 
be relaxed if the masses of the first- and second-generation squarks are large (> TeV) 
while the third-generation masses remain relatively small (< TeV), or in the presence of 
cancellations between the contributions of different CP-violating phases. We note that 
Ml also contributes at the one-loop level, but again, if accidental cancellations are allowed 
between terms, it remains essentially unconstrained. The phases of the third-generation 
trilinear couplings, (pAt b r have weaker constraints, as they contribute to EDMs only at 
the two-loop level. Again, accidental cancellations can occur that weaken further the 
constraints: see [7-17]. A comprehensive summary of the EDM constraints and other 
CP- violating effects in SUSY is given in [18]. Here we study the complete range of CP 
phases in order to sec the general dependences exhibited by our observables, and what 
luminosity might be required to observe these within the LHC environment. Therefore, 
we do not calculate explicitly which values of the other phases might be required for the 
points in our displayed scenarios to satisfy the EDM constraints. 

The precise determination of these phases is expected to be possible only at an e+e~ 
linear collider, for instance at the planned International Linear Collider (ILC) or at the 
Compact Linear Collider (CLIC). However, it will be crucial for such future search strate- 
gies to use LHC data to learn as much as possible, as early as possible. Furthermore, 
the combination of independent measurements at the LHC and a linear collider will be 
important to determine the Tinderlying model. 

In this paper we concentrate on the potential for observing unique CP-violating effects 
in decay chains at the LHC, and investigate the circumstances under which a determina- 
tion of the complex MSSM phases may be achievable, possibly with the support of other 
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LHC measurements. 

Specifically, we consider the LHC process gg — > it, with subsequent decay i — tx2, 
X2 ^"""^"Xi- We consider the situation where the X2 decay is a three-body decay; this 
leads to CP violation as there is a non-negligible contribution from interference diagrams. 
This process involves the three phases 0Mi, 0/^ and (pAt', we discuss below the combinations 
to which this process is sensitive. We extract information on the phases using triple 
products formed from the decay products of the stop. Such T-odd variables have also 
been studied in the context of heavy squark and stau decays in [6,19-23]. Other related 
studies are [24-34]. 

The first CP-odd asymmetry we consider is formed from 7^ = pt ■ {pi+ x p^-). This 
quantity has been studied at the parton level in [6], assuming pure gaugino-like neutrali- 
nos. In our current study we provide analytic expressions for the squared amplitude of 
the cascade process including full spin correlations and general neutralino mixing, and 
also provide an analytic expression for the phase space in the laboratory system. We also 
incorporate parton density functions (pdfs) and discuss the CP-odd observables at both 
the parton and the hadronic levels. Transition to the latter level has a big dilution effect 
on the measurability of a CP-odd asymmetry. We include the possible LHC uncertainties 
in masses and asymmetries and discuss the extent to which CP- violating phases may be 
constrained in such cascade decays at the LHC. 

In [21], further CP sensitive asymmetries formed from the momentum of the b quark 
in the top decay were studied under the assumption of 2-body neutralino decays into 
on-shell sleptons, namely %, = py, ■ (p^+ x p^-) and %b = Ph ■ (pt x P£±)- These variables 
are sensitive to and (pAt, but have different dependences on the CP-violating phases 
as described in Section 12.31 Therefore, a combination of all three observables would in 
principle allow one to disentangle the influences of all three phases. 

Since T-odd observables can also be generated by final-state interactions at the one- 
loop level, one should in principle combine the asymmetry for a process with that one of its 
charge-conjugated process. If a non-zero asymmetry is then observed in this combination, 
it must correspond to a violation of CP symmetry. For the triple product, 7^, this is 
experimentally possible as long as the associated W decays into a final-state lepton, 
which enables us to determine the change of the t. Regarding the other triple products, 
we require information from the opposite decay chain to identify the charge. In all the 
scenarios we consider, the decay t xt^i dominant, enabling charge identification in 
principle. However, a detailed simulation including all combinatorial aspects and also 
other background processes would be required to validate this possibility. 

We begin by describing the process under consideration in Section [21 including the 
phases involved and their various effects. In Section [3] we present numerical results for 
three specific benchmark scenarios and discuss the potential for a measurement at the 
LHC. The Appendices contain details of the Lagrangian, the expression for the squared 
amplitude including full spin correlations, and the kinematics of the phase space in the 
laboratory system. 
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2 Formalism 



2.1 The process studied and its squared amplitude 

We study the dominant stop production process at the LHC, namely 

gg iiii, (1) 

with the subsequent decay chain 

U^x',+t^xT(^' + Wb. (2) 

At tree level, the production process ([1]) proceeds via g exchange in the direct channel and 
t exchange in the crossed channel, and via a quartic coupling, as shown in Fig. [H Another 
possible source of tis is their production in gluino decays, g tt. However this leads to 
an experimentally more complex topology than the direct production and consequently 
we do not investigate this channel. The Lagrangian and the resulting neutralino and stop 
mixings and couplings are described in Appendix \^ 

Since gluons do not couple to off-diagonal combinations of stop mass eigenstates of 
opposite chirality, and similarly for stop exchange and the quadratic couplings, tit2 pro- 
duction occurs only at the loop level, and we do not consider it here. We focus here 
on till production, since the reconstruction of full decay chains of ti seems achievable, 
even in the complex experimental environment at the LHC. With the exception of the stop 
mass eigenvalues, see Appendix [Aj no effects from supersymmetric CP- violating couplings 
occur in the tree-level production process. 

The first step in the cascade decay chain is the two-body process U tx2- Here 
CP- violating couplings of the ii enter as well as those of the X2; ^^^1 are dominated 
by the phases (pAt and 0Mi, see Appendix IB. II \ Since constraints from electric dipole 
measurements constrain strongly 0^, we set 0^ = in our study [35-37]. We consider 
spectra where the second steps in the cascade decay chains are the three-body decays of 
the neutralino, X2 Xi^~^^~ (cf- Appendix IB. 2|) and the dominant top decay t — > Wb (cf. 
Appendix IB.3p . The neutralino decay occurs via exchange in the direct channel and 
via ii^R exchanges in the crossed channels, cf. Fig. [2J It is very sensitive to CP- violating 
supersymmetric couplings, and its structure has been studied in detail in [38,39]. The 
phase (pMi (and also 0^, which has been set to zero here) affects the mass of the X2> 
well as its couplings and decay rates. 



Using the formalism of [39,40], the squared amplitude |Tp of the full process can be 
factorized into the processes of production gg — > titi and the subsequent decays ii tx2, 
X2 Xi^^^~ and t Wb, with the second ii being unobserved. We apply the narrow- 
width approximation for the masses of the intermediate particles, ti, X2 ^^^d t, which is 
appropriate since the widths of the respective particles are in all cases much smaller than 
their masses, cf. Table 13. 1[ the squared amplitude can then be expressed in the form 

3 

\T\' = 4\A{h)\'\A{xl)\'\A{t)\'P{hh){P{^^^^^ 

a=l 

^ Their structure has also been studied in detail in [21]. 
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Figure 1: Feynman diagrams for the production process gg titi. 




Figure 2: Feynman diagrams for the three-body decays x° Xk^~^^ ■ 



+ E ^pit)^D{t)D{xl) + E ^pixlt)^D{xl)^D{t)h (3) 

6=1 a,b=l 

where a = 1,2,3 refers to the polarisation states of the neutrahno Xi* and top quark, 
which are described by the polarisation vectors s"'{x^), s^{t) given in Appendix IB.li In 
addition, 

• A(ti), A(x2) and A(t) are the 'propagators' of the intermediate particles which lead 
to the factors Ei^/miJ^i^, i^^o/mj^oF^o and Et/rritTt in the narrow-width approxi- 
mation. 

• P{iiii), P{tx2), D{Xi) and D{t) are the terms in the production and decay that are 
independent of the polarisations of the decaying neutralino and top, whereas 
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• Sp(x°), Sp(t), Ep (x2^) and S^(xi'), ^^^(t) are the terms containing the correlations 
between production and decay spins of the X2 t. 

According to our choice of the polarisation vectors ['S^(t)], see eqs. ( I60|) - (!65l) in 

Appendix lB.il Sp/P is the longitudinal polarisation, Sp/P is the transverse polarisation 
in the production plane, and Sp /P is the polarisation perpendicular to the reference plane 
of the neutralino Xi [^op quark t]. 

2.2 Cross section for the whole process at parton level 

The differential cross section in the laboratory system is 

dcr = ^|T|2(27r)V(pi +P2 - 5^p„)rflips(p„), (4) 

^ n 

where is the beam energy of the gluons, pi and p2 are the momenta of the incoming glu- 
ons, the p„ are the momenta of the outgoing particles and d\ips{pn) is the Lorentz-invariant 
phase space element. Integrating over all angles, all spin-dependent contributions are can- 
celled and the cross section for the combined process of production and decay is given 
by 

a = a{gg iJi) x BR{ii txl) x BR{xl ^ Xi^^i') x BR{t Wh) 
^ |A(tOp|A(x^)nA(t)p 
2EI 

X I P{t,h)P{xlt)D{xl)D{t){2TiY5\p,+p, - Y,Pn)dYiV<Pn). (5) 

n 

The explicit expression for the phase space in the laboratory system is given in Ap- 
pendix o 



2.3 Structure of the T-odd asymmetries 

Suitable tools to study CP-violating effects are T-odd observables based on triple products 
of momenta or spin vectors of the involved particles. In this paper we study the following 
T-odd observables: 

% = Pf {Pe+ X pi-) , (6) 
% = Pb- iPe+ X , (7) 
%b = Pf {Pb X pt±). (8) 

The T-odd asymmetries are defined as 

A _ iVr,+ -iVr,- _ /sign{r,}|TPrflips . _ , . tfe f9l 

where Nr^+i ^Tf- ^'^^ the numbers of events for which Tf is positive and negative respec- 
tively, and the second denominator in ([H]), /|Tp(ilips, is proportional to the corresponding 
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cross section, namely cr{gg — > titi tx^i^i^) in eq.([6]) and cr{gg iiii — > Whx^d-^ ) 
in eqs. (I7j) and ([H])- In the second numerator in ([H]), only the triple-product correlations 
enter via the spin-dependent terms, as explained in eq. (ITT]) and the following sections. 

The observable At^ has the advantage that it is not necessary to reconstruct the mo- 
mentum of the decaying t quark. However, as explained below, in order to disentangle the 
effects of both phases of At and Mi, it will be necessary to study all possible observables. 

The asymmetry Atp eq. ([9]), is odd under the naive time-reversal operation. It is the 
difference of the number of events with the final top quark or 6-jet above and below the 
plane spanned by x pi- in eqs. ^ and ([7]) , and by x f)i± in eq. ([H]) , normalised by 
the sum of these events. 

As can be seen from the numerator of Atj , in order to identify the T-odd contributions, 
we have to identify those terms in |Tp, eq. ([3]), which contain a triple product of the form 
shown in eqs. ([SD^dH])- Triple products follow from expressions ie^^p^a'^b'^c^'d'^ , where a, 
b, c, d are 4-momenta and spins of the particles involved, which are non-zero only when 
the momenta are linearly independent. The expressions iep^ypaa^b^c^d'^ are imaginary and 
when multiplied by the imaginary parts of the respective couplings they yield terms that 
contribute to the numerator of Axf, eq. ([9]). In our process, T-odd terms with e-tensors 
are only contained in the spin-dependent contributions to the production, Sp and in 
the spin-dependent terms in neutralino decay, S^(x°). It is therefore convenient to split 
Hp (x°t) and into T-odd terms Sp'°(x°t) and S^*^(x°) containing the respective 

triple products, and T-even terms Sp'^(x°t) and without triple products: 

^p{x]t) = s^^'°(x?t) + s« (x?) = + ss'^lx?). (10) 

The other spin-dependent contributions Sp(xj) and Sp(t), as well as S|)(t), are T-even. 

When multiplying these terms together and composing a T-odd quantity, the only 
terms of |Tp, eq. ([3]), which contribute to the numerator of Axf are therefore 
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The first term in eq. (flT]) is sensitive to the T-odd contributions from the production of 
the top and the neutralinos x'j- Comparing eq.([HZD with eqn. (17^ and (1^ leads to the 
following possible combination of contributing momenta 

S^^'0(x°t)S^^(x?)S5,(t) ~ 6,.,.."'^(x;)p^oa^'^(t)K X {P[i.,e-]sn{Plt,w]s'). (12) 

The second term and third terms in eq. fill I) are only sensitive to T-odd contributions 
from the neutralino x'j decay. The second term depends only on the polarization of 
Xj, comparing eq. fl86|) with eq. fl57|) therefore leads to the only possible combination of 
momenta 

Ef{x'j)^f{x',) ~ {pts^ X 6,.,.s>^op^p^+. (13) 

Since the third term depends on the polarization of both fermions, Xj and t, the possible 
combinations, comparing eq. fl86l) with eqn. fl66l) and fl92l) . are 

E^^'^(x;t)SS°(x;)S5,(t) ~ (p*s'')(Pxr')(^V[.,iy])xe^.p.s>^„p^p^+ (14) 

and {s''s''){s''p[t,w]) x e^,p,s>^op^_p^+. (15) 
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As can be seen by substituting eqs. dSOD-dSSD into eq. ([67D in Appendix EH Ep'^(x°t) 



vanishes for the combinations (ab) = (11), (22), (33), (13), (31), because they contain cross 
products of three hnearly-dependent vectors. Only for the remaining combinations, (ab) = 
(12), (21), (23), (32) , do we get a T-odd contribution to the production density matrix. 

Similarly, the expression for the T-even contributions, Sp''^(x^t), eq. fl66l) in Appendix 
IB.lt has non-zero components for a = 1,3 but vanishes when a = 2. These expressions 
are multiplied by S^°(x°), eq. fl86l) . and therefore only S^°(x°) and E^'^(xj) contribute. 

In the following section we derive the three triple products, study their different de- 
pendence on phases and provide explicitly a strategy for determining (pAt and (pMi and 
disentangling their effects. 

2.4 Strategy for determining (f)At and 
2.4.1 Derivation of the triple products 

In order to describe the spin of a fermion / in general, we introduce three four vectors, 
s^(/), a = 1,2, 3, such that the s"" and the momentum and mass of the fermion p/m form 
an orthonormal set of four- vectors [40]; 

p-s" = 0, (16) 
s^.s^ = (17) 

s^^st = -g,. + (18) 

where repeated indices are implicitly summed over. 

Applying eq.(|T8l) on eqs.f|T2l)-(|T5l) lead to kinematic expressions that contain only 
explicit momenta. Expanding terms with tf^^pa in time- and space- components gives 
scalar triple products between three momenta. 

In our process we can classify the terms of eq. lITT]) as follows: 

• The terms of eq. (fT2|) lead to a combination between Ttb and Tf,. 

• The terms of eq. f|T3|) lead only to %. 



• The terms of eq. (fT^ lead again only to % but terms of eq. (fT5|) produce % as well 
as %, due to interference effects between both spin vectors of pt and p^o. 

2.4.2 T-odd terms sensitive to % 

We consider first %, eq. ([6]). As this includes the reconstructed top quark momentum, 
there are no spin terms from the decay of the top quark and the contributing terms are 
the second and third term in eq. fllip as explained in the previous paragraphe. The CP- 
sensitive terms of the decay density matrix are given by eqs. and the contributing 
kinematical factor is g2, eq. flQOl) . 



9t = ^mfce^,p,s">^op^p,% . (19) 
We note that is purely imaginary. When inserted, for instance, in eq. 0881) it is mul- 



tiplied by the factor i ■ Ira{fhf^^O'l^*}, which depends on the phases 0Afi (and 0^) and 
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contributes to . Analogous contributions follow from eqs. (l87Il and (!89l) . The cor- 
responding T-even terms of the production density matrix also entering in eq. (ITT]) are 
obtained from eq. (IFTIl . 

2.4.3 T-odd terms sensitive to Ttb 

For the triple product T^f,, eq. ([8]), only the first term in eq. ffTTl) contributes, but the 
kinematics is complicated by the fact that we need to include the decay of the t in addition 
to that of the X2- This comes from the fact that the kinematical term that generates the 
triple product is ff', eq. fl^ : 

ft = e,,p„s'^^>^{x'l)plf'''m- (20) 

As both and s^'''{t)p1 are contained in this term, we need to include their decays 

in order to to produce a non-zero contribution. 

This term occurs only once in the t decay amplitude, eq. fl671) . and is multiplied by 
the complex pre-factor g^Im{aijb*j), eq. fl69l) . Both a^j and bij contain terms from the i 
and Xj mixing matrices, and so are sensitive to both the phases (j)At and (and 0^). 

2.4.4 T-odd terms sensitive to % 

The triple product %, eq. (^^, is the most complicated, as it contains contributions from 
both the t and X2 decays (the first and third terms in eq. (ITT!) ). The kinematics is rendered 
more complex by the need to multiply each T-odd contribution by the terms from the 
other two decays. Each T-odd component is generated through g2 and ft, as for the 
other two triple products. As a consequence of having a dependence on both the i and 
X2 decays, % is also sensitive to both phases (pAt and (f)Mi (and 0^). 

2.4.5 Disentangling of effects of (pAt and (pMi 

The T-odd asymmetries, eq.(l9]), are determined by those CP- violating couplings that are 
multiplied with the respective triple product. Under the assumption that (p^ is small, 
the neutralino sector depends only on 0jv/i and the stop sector only on (pAf Since the 
involved triple product momenta show different dependence on the CP-violating phases, 
as discussed above, it is possible in principle to disentangle the effects of (pAt and (pMi in 
our process and to determine the phases separately. 

The decoupling is possible as the triple product %, can only be produced by the term, 
Sp^(Xj')S^*^(x°), cf. section [2.4.11 The T-odd contribution in this term comes from the 
decay of the Xj and consequently is only sensitive to the phase (pMi- Once we have used 
the triple product % to determine the phase (pMi we can then use the value as an input 
for the triple products, %h and %, in order to determine the phase (pAt- 
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3 Results 



3.1 Scenarios 

In this Section we analyse numerically the various triple-product asymmetries introduced 
in eqs. (I6])-® at both the parton level and with the inclusion of parton distribution 
functions (pdfs) to study the discovery potential at the LHC. In particular, we study 
the dependences of these triple-product asymmetries on the MSSM parameters Mi = 
|Mi|e*'^*'^i and At = \At\e^'^^-i- . We also analyse the effects of these parameters on the 
masses and branching ratios of the particles involved in our process. 

For our numerical analysis we study in detail at both the partonic and pdf level a 
reference scenario, A, where the x5 is a gaugino-higgsino mixture. For comparison, we 
also study at the partonic level a non-universal Higgs masses (NUHM) scenario, B, and a 
third scenario, C, in which the x\ is higgsino-like. The particle spectra for these scenarios 
have been computed with the program SPheno [44]. These three scenarios have been 
chosen to have similar masses, as displayed in Table 13.11 so that the kinematic effects are 
similar in each case. We perform our studies using our own program based on the analytic 
formulae we have derived for the various cross sections and spin correlations. The program 
uses the VEGAS [45,46] routine to perform the multi-dimensional phase-space integral. We 
constrain ourselves to cases where m^o < m^o + rrizo and m^o < rrin , so as to forbid the 
two-body decay of the Yj- The branching ratios for both processes have been calculated 
with Herwig++ [47,48]B 

The feasibility of measuring these asymmetries at the LHC depends heavily on the 
integrated lumijiosity at the LHC. For this reason we look closely at the cross section, 
a = a{gg iih) x BR{ii txl) x BR{x2 ^ X?^^^") x BR{t Wb) and determine 
the nominal luminosity required to observe a statistically significant result. 

3.2 CP asymmetry at the parton level 
3.2.1 Dependence of m^o and At on (pMi and (pAt 

We start by discussing the dependence on Mi = iMile*"^*^! of the parton-level asymmetries 
for each of the three scenarios. In order to see the maximum dependence upon 0Mi, we 
use the reconstructed t quark momentum and the triple product % = pt ■ {p£+ x p^-) . 
It should be noted from the following plots that the asymmetry is obviously a CP-odd 
quantity that in addition to a measurement of the phase, also gives the sign, as seen 
in Fig. [3](a). In comparison, using CP-even quantities, for example the mass, it is not 
possible to determine if the phase is positive or negative, as seen in Fig. [3](b). 

We see in Fig. [3](a) that the biggest asymmetry appears in scenario A, which attains 
|»4.Tt|max ~ 12% when (pMi ~ O.Stt. One aspect of the plot that may be surprising is that 
the asymmetry is not largest at the maximal value of the phase {(pMi = f )• This is due 
to the coupling combinations and interferences and can be seen from the equations in 

^Beyond the Standard Model physics was produced using the algorithm of [49] and, in the running of 
ttEM, the parametrization of [50] was used. 
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Scenario 


A: Reference 


B: NUHM - 7 


C: Higgsino 


Ml 


109 


97.6 


105 


M2 


240 


184 


400 


fi 


220 


316 


-190 


tan P 


10 


20 


20 


Ml 


298 


366 


298 


Me 


224 


341.7 


224 


Mq3 


511 


534.5 


511 


Mu3 


460 


450 


460 


At 


-610 


-451.4 


-610 


Mi, 


396.5 


447.8 


402.6 


Mi, 


595 


609.6 


591.6 




177 


172.8 


186.3 




301.6 


346.05 


421.1 




302.4 


369.8 


303.1 




229.2 


345.2 


229.2 


0Mi 





TT 





TT 





TT 


X-l 


100.8 


106.1 


94.8 


96.3 


99.2 


97.6 




177.0 


171.3 


167.1 


166.6 


186.2 


179.8 


X-'i 


227.9 


231.8 


323.8 


325.5 


199.4 


206.2 




299.1 


297.6 


343.4 


341.8 


419 


418.9 




3.88 


3.88 


3.48 


3.48 


5.29 


5.29 




1.4x10-^ 


1.4x10"^ 


2.3x10^5 


2.3x10"^ 


3.0x10-3 


3.0 xlO-3 



Table 1: Parameters and spectra for the three scenarios A, B, C considered in this paper. 
We display the input parameters \Mi\, M2, tan/3, m^^ and m^^ and the resulting 
masses rn^o, i = 1,...,4, for (^Afi,^^) = (0.57r, 0) and (0.57r, O.Svr). The parameters 
M2, |/i| and tan/? in scenario B are chosen as for the scenario SPSla in [41]. We used 
mt = 171.2 GeV [42] and the SM value for the top width ~ 1.5 GeV [43] for our study. 
All masses and widths are given in GeV. 



Sections 12. 4. 2[ In Fig. [3](b), the dependence of the masses of the neutralinos is shown. It 
can be seen clearly that the variations are too small to be used to determine the phase. 

In the cases of the two other scenarios shown in Fig. [3](a), the dependence of the 
asymmetry on the phase 0Mi is similar but slightly smaller. In the case of scenario B 
(NUHM), the peak asymmetry is I^Tjmax ~ 9% when (pMi ~ 0.37r and in scenario C 
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Figure 3: The asymmetry at threshold for the production process gg — > tt for scenarios 
A, B and C for (a) Art ^ function of 4>Mi, and (b) the masses of the neutrahnos as 
functions of (pMi- 
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Figure 4: (a) The asymmetry Art, threshold for the production process gg — > tt for 
scenarios A, B and C, and (b) the asymmetry A% at threshold, both as functions of (pAf 



(Higgsino) it is l^rjmax ~ 7% when (pM^ ~ 0.257r. Again, the asymmetry does not peak 
when the phase is maximal. 

To study the dependence upon (j)At we need to use the triple products sensitive to this 
phase, % = pf {pe,+ x p^- ) and Ttb = Pf {Pb x Pe± ) • Fig. Hl^a) shows Ar,, and we see that the 
biggest asymmetry again occurs in Scenario A, but the maximal asymmetry is only about 



half of |.4rt|max with \A' 
Scenario A, with {At^ 



i, I max 



rjmax ~ 6%. Scenario C produces a very similar asymmetry to 
5.5%, whereas the asymmetry in Scenario B is much smaller: 
l-^rjmax ~ 2.5%. Fig. 111(b) shows that the general shape of the asymmetries for Ar^^ is 
similar to that of Ar^ apart from a difference in sign and that all the asymmetries are 
actually slightly larger. In fact, for Scenario C, the largest asymmetry is generated using 



Ttb with Ar. ~ 8% when 



O.Svr. 



In the subsequent analysis, we concentrate on the favourable Scenario A, with just a 
few remarks on the others. 
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3.2.2 Contour Plots of At^ and At^^ for Variable Mi and At 
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Figure 5: Contours in scenario A (in %) of the parton-level asymmetries (a) Art for the 
triple product % = pt ■ (p£+ x p^-), as functions of the variables Mi and (pM^, and (b) 
Axti for the triple product Ttb = Pi+ ■ (pt x pb), as functions of the variables Mi and (pAf 

If we now lift the restriction of the GUT relation for |Mi|, we can see how the asym- 
metry varies with |Mi| while leaving all the other parameters the same, for Scenario 
A. Fig. shows that the asymmetry peaks at |Mi| ^ 130 GeV and ~ 0.257r 

when \ATt \ ~ 15%. Importantly though, the asymmetry can remain above 10% between 
|Mi| = 110 GeV and |Mi| = 190 GeV, which is most of the range allowed in this scenario. 

By including the decay of the t quark that was produced in the t decay, we can 
also study the effect of (pAt on our asymmetries. As the spin-correlation information is 
now carried by the t quark, we have to change the triple product used to measure the 
asymmetry, eq. (jH])- It is found that the largest asymmetry can be measured using the 
triple product, Ttb = Pe+ ■ (Pt xPb) where I^TtJmax ~ 8.5% when (pAt ~ 0.57r in Scenario A, 
as seen in Fig. [5](a). It may be noted that this asymmetry is slightly smaller than those 
of [21] that can be reconstructed experimentally. In that paper scenarios were chosen 
where the X2 decays via a two-body process, whereas here we concentrate on scenarios 
where the X2 decays via a three-body process, so to maximise the sensitivity to 0Mi- This 
phase dependence can also be seen with the triple product % = p^- (p^+ x p^+) although 
the asymmetry is found to be smaller here with |.4Tjmax ~ 6%, see Fig. [6](a). 

We have also considered the dependence of the asymmetry on a common trilinear 
coupling, A = At = Ab = A^, in scenario A, as shown in Fig. Mjo). It can be seen that 
the asymmetry is stable for the bulk of the region scanned, and only decreases near the 
edge of the acceptable region for our scenario. The peak is now |w4,TtJmax ~ 9%, when 
At ~ —500 GeV, and the region where |^Ttbl > 8% extends from At ~ —650 GeV to 
At ^ -250 GeV. 

We now consider the effect on the asymmetry of varying simultaneously both the 
phases 0Mi and (pAf The triple products % = pb- {pi+ x -p^-) and Ttb = Pf {pb x P£±) can 
have contributions from both phases, so we concentrate on these. For Fig. [Tl^a) shows 
that the area of parameter space where 0Mi and (pAt constructively interfere is actually 
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Figure 6: Contours in scenario A (in %) of the asymmetries (a) Ab for the triple product 
T = pb ■ (p^- X p^+), as functions of the variables Mi and (pAt and (b) Atb for the triple 
product T = p^+ ■ (pt X Pb), as functions of the common variables A = At = Ab = Ar 
and (pAf 



Ar^, Scenario A, ^/I ~ 2mj 



a 



0.5 



1 



1.5 



1 




\ 

\ 


1 / 

t .' .' 

! ' .' 




i' 

-2/ 


i 

/ 


\ \ • 

■\. 






0/ 


'a 


■ ■. \ 

4'. i d 
; 1 
/ / 




1 

i 
i 
i 
i 


/2 /' /' ^■ 

.'4 /' 
, I 6 


.' / 

■' .-■ / ' " 
/ /' .'' 
/' ,■' 


-2'\ 


\ 
\ 
\ 


I 



0.5 



1.5 



Axt^, Scenario A, ~ 2mj 
(b) , , r- 



0.5 



t/TT 1 



>8 ^ H 12 '. 1 
j ■ 



1 ' 2 ! 4 ; 6 ; 



1 .5 *8 i-eU ';-2 ) -1 



2 t.^^ 1_ 

0.5 




-1 -2 ■, -4 1-6 1,-8: 



1 1.5 2 

0Mi/71" 



Figure 7: Contours (in %) of the asymmetry at the parton level in scenario A with 
Ml = 109 GeV for the triple products (a) % = Pb- (P£+ x p^-) and (b) T^b = p^-i- ■ (pb x pt) 
for varying phases 0Mi and (pAf 



quite small and peaked around 0Mi ~ 0.27r and (pAt ~ O.Svr. Apart from this area, varying 
both phases generally results in a reduction in the asymmetry observed, caused by the 
neutralino and squark mixing entering the couplings. Section I2.4.3[ Importantly when 
(pMi ~ TT or (j)At ~ the asymmetry vanishes, as it should. Fig. [7](b) demonstrates that, 
for this scenario, (pMi generates virtually no asymmetry for %. However, (pMi can still 
significantly reduce the asymmetry that cpAt can produce and, again, when cp^i ~ tt we 
see that {Axt^l ~ as expected. 

If we now modify scenario A slightly by setting |Mi| = 160 GeV, this results in a more 
interesting scenario as the phases 0Mi and (pAt can interfere constructively to produce an 
asymmetry larger than that seen before. When 0jv/i ~ 0.47r and (pAt ~ l.Svr, we observe 
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A%, Scenario A, \/3 ~ 2mj 
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Figure 8: Contours (in %) of the asymmetries in scenario A at the parton level with 
Ml = 160 GeV for the triple products (a) Tf, = Ph {Pe+ xp^-) and (b) Ttb = P^+-(PbXPt), 
as fucntions of the varying phases 0Afi and (pAf 



a peak asymmetry, \AtJ ~ 7% for the triple product Tf,, as seen in Fig. [HI 



3.3 Dependences of branching ratios on cpM^ and 



BR{ii ^ xlt), Scenario A (0^, = 0) BR{xl x\i^t~), Scenario A (0^, = 0) 
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Figure 9: Contours (in %) of branching ratios in scenario A as functions of Mi and (pMi'- 
(a) BRik ^ X^t) and (b) i?i?(x^ ^ xTH, ^ = e or /i. 



In order to determine whether an asymmetry could be observed at the LHC, we 
need to calculate the cross section for the total process. Important factors in the total 
cross section are the branching ratios BR{ti —>■ x^^) (fo^ CP- violating case see [51]) and 
BR{x2 ~^ [38]. Both of these change considerably with and (pAt, altering 

the statistical significance of any measurement of I^t^I- Analyzing first the variation 
with Ml, seen in Figs. 19) and fTOlfa). we see that the branching ratio BR(ti — >• X2^) is 
indeed sensitive to variation of the phase, but can vary more strongly with |Mi|. For 
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Figure 10: Contours (in %) of the branching ratio BR(ti X^'^), in scenario A as 
functions of varying (a) Mi and (pAt and (b) the common trihnear couphng At = = A^ 
and the phase of the top-quark trihnear couphng (pAf 



example, if (pMi = tt when \Mi\ ^ 150 GeV then BR ^11%, but if we keep the phase 
the same and change to \Mi\ ^ 180 GeV then BR ^ 4% (i.e., it drops by almost a factor 
of four), as seen in Fig. [9](a). The general reduction of BR{ti X2't) as Mi increases is 
to be expected as the character of X2 '^iH be less gaugino-like. Similar large differences 
are found in BR{j(^ — >• which varies between 3% for Mi < 135GeV and 9% for 

Ml ^ 165GeV Fig. Mh). 

The phase 0At does not enter BR{x2 ~^ but can have a large effect on 

BR{ti — i> X2^)- III scenario A, we see in Fig. [10] that BR ^ 8% at (pAt = but increases 
to BR ^ 24% at (pAt = ^ (i-e- a factor of 3 increase). The branching ratio BR{ti — > Xq,^) 
also has a dependence on \At\ an this is shown in Fig. [TU](b). We see that if (pAt = 
then the branching vary between, BR ^ 4% when \ At \ ^ — 750GeV and BR ^ 12% when 
\At\ ^ -lOOGeV. 

In the range of Mi = \Mi\e''l"''i and At = \ At\e''*'^t studied, we find that BR{ti xl^) 
varies between 4% and 24% and BR{x2 ~^ Xi^^^~) between 2.5% and 9% for scenario A. 
Similar plots can also be produced for scenarios B and C but are not presented here. It 
is found that BRiJi xV^) varies between 4% and 14% for scenario B and between 8% 
and 35% for scenario C. For BR^y^ Xi^"*"^") the variation is between 3% and 12% for 
scenario B and between 2% and 5% for scenario C. 

3.4 Influence of parton distribution functions (pdfs) on CP asym- 
metries 

So far we have studied the triple-product asymmetries only when the production process 
is close to threshold, and the ti pair is produced almost at rest in its centre-of-mass frame; 
triple-product effects due to spin effects are usually greatest close to threshold. However, 
production at the LHC is not in general close to threshold, and we must include pdfs in 
our analysis to see how an initial boost to the ti affects the asymmetry. We focus on 
scenario A: similar results are obtained in Scenarios B and C. 
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Figure 11: (a) Asymmetry for scenario A as a function of the t momentum, (b) Total 
cross section for scenario A for gg — > tt as a function of the parton-parton centre-of-mass 
energy. 

Fig. fTDfa) shows the asymmetry I-AtJ as a function of the ti momentum, and shows 
clearly that the asymmetry is peaked at the threshold for ti production, where the stops 
are produced almost at rest, and that it falls sharply as the energy increases. Fig. [TlTb) 
shows the total cross section in 14 TeV collisions at the LHC for gg — »• titi as a function 
of the parton-parton centre-of-mass energy, and demonstrates that the peak production 
occurs close to threshold with a long tail of production at high energy. In addition, even 
when production occurs at a low parton-parton centre-of-mass energy, in the majority 
of cases one gluon may be carrying significantly more momentum than the other in the 
collision. Consequently the produced ti can have a large longitudinal component to its 
momenta. Both these factors mean that the asymmetry observed at the LHC will be 
substantially smaller than if the all ii were produced at threshold B In should be noted 
that similar results were found for all asymmetries and scenarios, and this 'dilution' factor 
is always present. 

We use the MRST 2004L0 pdf set [52] in our analysis of the asymmetry, and plot the 
integrated asymmetry l^rj as a function of 0Mi and at the LHC in Fig. [T^ a). as the 
solid line. We see that the inclusion of the pdfs reduces the asymmetry by about a factor 
of four in this case. This reduction is unsurprising, given the reduction in asymmetry 
when one moves away from threshold shown in Fig. [TlTa). though the dilution factor does 
depend on the scenaric0. 

Using the production cross sections and branching ratios we can then estimate the 
integrated luminosity required to observe an asymmetry at the LHC. We assume that 
Nr±, the numbers of events where T is positive and negative as in eq. (Q, are binomially 

^Both these effects could be overcome if one could measure the stop-stop invariant mass and tag the 
stop momenta, but this is unlikely to be possible with great accuracy. 

■^These results have been checked independently using Herwig++ [47, 48] with three-body spin corre- 
lations included, a feature that is currently not available in an official release of the code, but will be 
included in a future version. 
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Figure 12: Integrated asymmetries with parton density functions included in the pro- 
duction process. The dotted and dashed hues indicate the asymmetry required in order 
to observe a la deviation from zero with the indicated luminosities, see the text: (a) 
% = Pt ■ {pe+ X P^-) in scenario A as a function of (pMi with Mi = 130 GeV, and (b) 
%b = Pe+ ■ (Pt X Pb) in scenario A as a function of (pAt with Mi = 109 GeV. 




Figure 13: Integrated asymmetries with parton density functions included in the pro- 
duction process. The dotted and dashed lines indicate the asymmetry required in order 
to observe a la deviation from zero with the indicated luminosities, see the text: (a) 
Tfc = Pb ■ {p£+ X Pi-) in scenario A as a function of (pAt with Mi = 109 GeV. 



distributed, giving the following statistical error [53]: 

A{AtT^' = 2^/eil-e)/N, (21) 

where e = Nr+ / {Nr+ + Nt-) = \{1 + At), and N is the number of selected events. This 
can be rearranged to give the required number of events for a desired significance. 

Figs. [12] (a), (b) and fT^ a) show the expected levels of the integrated asymmetries in 
scenario A with pdf effects included (solid line) together with dotted and dashed lines 
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showing the level of asymmetry one would need with the corresponding integrated lumi- 
nosity in order to obtain a statistical error At > ^{At)- In other words, the asymmetry 
could only be seen at the level of la where the solid line is above the relevant dotted or 
dashed line. For example, in scenario A after 100 fb~^, the asymmetry could only be seen 
for a small area of parameter space around 0Mi = O.SSvr and I.Ttt. Figs. [12] (a) and (b) 
show that even if (pMi or (pAt has a value that produces a maximal asymmetry, we require 
a substantial integrated luminosity if we are to find a statistically significant result. In 
addition, it must be noted that we have not included any detector effects into our analysis, 
and one could expect that the required integrated luminosity would rise substantially after 
the inclusion of backgrounds, trigger efficiencies, etc. A measurement of an asymmetry 
with an accuracy of a few % might be possible with 100 fb~^ of integrated luminosity, but 
it would probably be insufficient to constrain significantly the model parameter space. 
However, an interesting measurement could be made with an integrated luminosity above 
300 fb~^, which is targeted by the proposed LHC luminosity upgrade. 

3.5 Determination of the CP-violating phases 

As we have shown, it will be challenging to determine the phases (f)Mi and (j)At in our 
process using the triple-product asymmetries alone. However, it would be very worthwhile, 
as a non-zero measurement of a T-odd asymmetry would provide unique evidence of CP 
violation. In the rest of this section, we examine briefly the potential for a measurement 
using other variables, again concentrating on Scenario A. 

3.5.1 Observables: masses, cross sections and CP asymmetries 
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Figure 14: (a) The mass of the stop squarks tj, j=l,2 as functions of 4>At/'^, and (b) 
contour plot showing the areas of the (Mi,0Mi) parameter plane consistent with a mass 
difference between X2 and Xi of 20, 40 and 60 GeV respectively. The bands assume a 1% 
error in experimental measurement of the mass difference and a 5% error in M2. 

Fig. m^b) showed how the masses of the x'^s vary with and Fig. [TWa) shows how 
the masses of the vary with (pAt in scenario A. Unfortunately, the variations in both 
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of these observables are only about 1 — 2 GeV, which are significantly smaller than the 
experimental errors expected for these measurements. A far more accurate measurement 
at the LHC will be the mass difference between X2 X15 fhis can be determined 
in our scenario with a clear dilepton end-point. The accuracy of this measurement is 
expected to be < 1% and, if we assume that M2 can be determined to 5% [54], we find 
the regions plotted in Fig. [TW b). At the smaller values allowed for Mi in scenario A, we 
see that this observable does not depend sensitively enough on (pMi for a measurement 
to become possible. However, as Mi increases we see that the sensitivity to (pMi becomes 
much clearer. Importantly, in scenario A, it is only possible to have a mass difference, 
X2 ~ Xi ^ 40 GeV if 0Mi is present. 

3.5.2 Inclusion of branching ratios 

Other observables sensitive to the phases 0Mi and (j)At are the branching ratios BR{x2 ~^ 
Xi^~^£~) and BR{ti — > X2^), as discussed in Section l373l As is the case for the masses, 
though, our current expectation of the accuracy of this measurement at the LHC looks 
insufficient to constrain the phases. Fig. [T^ a) shows in the context of scenario A that, if 
a measurement i?i?(x2 ~^ Xi^~^^~) = 0.4 is made and we assume that the accuracy at the 
LHC is 50% (Ai), then the constraints on Mi and are rather weak. However, if the 
accuracy could be improved to 10% (A2), a determination of Mi and (pMi looks possible 
if this analysis is combined with information from the Xi mass difference and that of 
the triple-product correlations. For the branching ratio, BR{ti —>■ X2^)) the conclusion 
is similar, as seen in Fig. [T5](b). With a measurement at 50% (Ai), we again see that a 
determination of the CP- violating parameter is not possible but, if a measurement can be 
made with an accuracy of 10% (A2), then a determination of (j)At would be more plausible. 



BRix^ ^ = 0.04, BR{h ^ m = 0.1, 

Scenario A {(j)At = 0) Scenario A (0Mi = 0) 




100 120 140 160 180 200 100 120 140 160 180 200 
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Figure 15: Parameter space allowed when the experimental accuracy of the branching 
ratio measurement is 50% (Ai) or 10% (A2) for (a) BR{x2 ^ = 0.04 and (b) 

BRik ^ x^t) = 0.1. 

Thus, we may be able to to pin the model parameters down with greater accuracy by 
combining information on the CP- violating asymmetries with this and other information. 
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4 Conclusions 



We have studied direct stop production followed by the decay ii t^^j X2 ~^ Xi? 
where the latter is a three-body decay and provide compact analytical expressions for 
the amplitude and phase space. We have specifically concentrated on measuring the 
CP- violating phases of the parameters Mi and Af. 

We have provided a thorough analysis of the contributions to this process which lead 
to non-zero asymmetries in the parameters %, % and %{, formed from triple products 
of reconstructable final-state particles. These are sensitive to different combinations of 
the phases mentioned above. We studied three spectra which had different neutralino 
characteristics at the parton level and also studied the (large) effect of including pdfs 
which had previously only been roughly estimated in the literature. 

We found that with the design integrated luminosity of the LHC of 100 fb^^, the 
statistical errors would probably remain too great to distinguish a non-zero asymmetry 
measurement from zero for most of the ranges of and (f)At, and we recall that this 
initial study did not include detector or background effects. However, with a luminosity 
upgrade, the accuracy will improve and it could be possible either to measure a non-zero 
value or else to provide limits on the possible phases. 

Triple products are not the only variables sensitive to the phases of the parameters. 
We found that a good measurement of the mass difference between the X2 Xi ^^u- 
tralinos could constrain significantly the (Mi,0Mi) parameter space. It is possible that 
measurements of the two branching ratios BR(ti — >■ X2^) and BR{x2 ^ Xi-^"*"-^") could 
also constrain both (f)Mi and (f)j^^, although this is heavily dependent on the experimental 
accuracy achieved. However, the disadvantage of both mass differences and branching 
ratios is that a non-zero value can potentially be faked by other values of the real param- 
eters. This is in contrast to the asymmetries from triple products which are uniquely due 
to CP violation. Therefore, even though these will be challenging measurements at the 
LHC, they are worthwhile experimental objectives. 
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Appendices 

A Lagrangian and couplings 

The interaction Lagrangian terms for the production processes are 

= g^d'^GX^UcGlG';,, (22) 
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The interaction Lagrangian terms for X2 decay are 



i^z°e,+i 
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cos Qw 



XmXn 2 COS Qw 

where the couphngs a^^ and 6^^ are given by 

n=l n=l 

Here TZl^ is the mixing matrix of the squarks and 
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Le = T^e - ei sin^ 9^^, Ri = -e^ sin^ G^f, 
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where Pl,r = |(1 T Ts); ^ = mt/{\/2m,w sm(3). Here, g is the weak couphng constant 
{g = e/ cos ©VK, e > 0), and T^i are the charge (in units of e) and the third component 
of the weak isospin of the fcrmion i, Qyy is the weak mixing angle and tan/? = V2/V1 is 
the ratio of the vacuum expectation values of the Higgs fields. The unitary (4 x 4) matrix 
Nmk that diagonalises the complex symmetric neutralino mass matrix is given in the basis 
(7,^0,^?,^,°) by [55]: 
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f Mie"^^icl, + M2sl, {M^e'^^h - M2)swcw \ 

{M^e'^^^ - M2)swcw Mie''t'^ic^ + M2S^ Mz 
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(37) 
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where the abbreviations sw — sm6w, cw — cos6w, S2/3 — sin 2/3, 02/3 — cos 2/3 have been 
used. 



The interaction Lagrangian for t decay is 

A.xg = 9Qi4kPR + biPL)xUi + h.c, (38) 
where the couphngs and 6^^. are given by 



2 2 

n=l n=l 



aj. = > ;(7ei)*^L, 4 = E(^i)*^L- (39) 



Here 7?.^^ is the mixing matrix of the squarks and 
with 

flk = -V2e,sin^H/iVfci- v^(T3,-e,sin2^,^)J^, (41) 

fkk = -V2e,sin^i^(tan^^iV*2-iV*i), (42) 

^ifc = Yt {Nl, sin /? - Nl, cos /?) (43) 

= {h^RkY- (44) 

We also use the following relations from [21]. The left-right mixing of the stop squarks is 
described by a hermitian 2x2 mass matrix which reads as follows in the basis {tL,tii): 



= -(^1, t\) I I I h (45) 

e'MM? I 

= M? + (i - ^ sin^ Qw) cos 2/3 m| + ml (46) 

= + ^ sin^ cos 2/3 m| + (47) 

= (^D* = ^*(^* - cot/3), (48) 
0t = arg[^t-//*cot/3]. (49) 

Here tan/3 = ^2/^1 with ^1(^2) being the vacuum expectation value of the Higgs field 
if°(i?2), TTT't is the mass of the top quark and Qw is the weak mixing angle, /i is the 
Higgs-higgsino mass parameter and Mg, Mfj,At are the soft SUSY-breaking parameters 

of the stop squark system. The mass eigenstates ti are (^1,^2) = {tL:tR)TZ^ with 
e*"^* cos 9i sin 6^ 

^* = I I , (50) 



M? 

tLL 




where 



— sin^f e "f'*cos^f 
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where 



-\Ml I Mj -ml 

cosOi = sinOt = (51) 

J\M? |2 + (m? -M? )2 ,/|M^ P + (m? - )' 

y ' tLR' ^ ti III' y -^-^ ^ 

The mass eigenvalues are 

m? =-f(M^ )^,/(M| -M| )2 + 4|M| pV (52) 

We note that we have (pi ~ (pAt fo^' l^tl ^ l/^l cot /3. 

B Explicit expressions for the squared amplitude 
B.l Neutralino production U x^t 

Here we give the analytic expression for the production density matrix: 

\MiU - = nx]t) + S?,(x;) + S^(t) + T.f{x]t), (53) 

whose spin-independent contribution reads 

2 

P(x]t) = |- {(|a,,f + |%r)(p*p^o) - 2m,m^oi?e(a.,6^)} , (54) 

where pt and p^o denote the four-momenta of the t-quark and the neutralino Xfc- The 
coupling constants can be simplified and shown to be 

(|a,,f + = |/.^/ + cos^ ei\fj + sin^ ^d/A/ 

+2sin%cos^,-/2e [e-'^-^h%{f^^ + ^O] , (55) 
i?e(a,,6^) = cos^ e-,Re{f^^h'^^) + sin^ 0,-i?e(/|;./i^^.) 

+i sin2^,-i?e [e'^nh%\' + e-^<^7i,/ft] ' (56) 



where is given in eq. (H9ll . 

The spin- dependent contributions are T-even and are given by 

Sp(x^) = ^ {(|%r - |a.,r)m,o(p,.'^(x°))} , (57) 



S^(^) = y |(l%r - |a.,nm,(p,o.''(t))| , (58) 

where s"'{x^) (s*(t)) denote the spin-basis vectors of the neutralino Xj (t-quark). Again 
the coupling constant can be simplified as 

i\hj\' - |a.,f ) = cos2ei\h%\' - cos^^d/i/ + sin^^d4r 

-2 sin cos %i?e [e-^<^*-/z*4. (4- + /l,-)] • (59) 
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The three spin-basis four- vectors s^, and form a right-handed system and provide, 
together with the momentum, an orthogonal basis system. They are chosen as: 



s iXj) = 0, ^ , (60) 

o'ix'') = 0, g^-^ , (61) 




The spin-system for the top quark has been chosen analogously: 

\ l(]9txp^o) xptiy 

s'{t) = 0, ,^ , (64) 

sHt) = -L(\p,l^pX (65) 
"^t V \Pt\ J 

and and E^q denote the energies of the top quark and the neutrahno x°, respectively. 

The terms that depend simultaneously on the spin of the top quark and of the neu- 
tralino can be split into T-even, Ep'^(x2^), and T-odd, Sp'*^(x2^), contributions: 

^P'^m = ^{2i?e(a,,6^)[(s'^(x?)p*)(5'(t)Px°) " iPtP^^)i^''ix'j)s\m 

+ mtm^o{s'^{x'j)s\m\a,j\' + \hj\')], (66) 
^f'^) = -g'lmia,jb*^)ff, (67) 
where the T-odd kincmatical factor is given by 

ft = e,.,,s^'^{x'j)p'-os''''^m, (68) 
and the coupling constant by 

Im{aijb*j) = cos' eilmU^h'^^) + sin' eilm{fj^^h%) 

+lsm29ilm [e^^*>*,/ + e-^^/i./^] • (69) 

B.2 Neutralino three-body decay x^- x\^^^~ 

Here we give the analytical expressions for the different contributions to the decay den- 
sity matrix for the three-body decay, where we sum over the spins of the final-state 
particles [39] . The contributions independent of the polarisation of the neutralino x° 

D{x]) = D{ZZ) + DiZh) + DiZin) + ^(44) + DildR), (70) 
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are given by 
D{ZZ) 

D{Zh) 



9^ 



cos^e 



w 



-\HZ)\\Li + R 



\0,'^\\g,+g2) + {ReO 



"L\2 
kj ) 



cos^ Qw 



LfRe 



L* A t* 

ik 



A'*ih)i20,fg, + 0,f*gs) 



+ /,^7.iA-(£i)(20;'f + Olfg,) 



D{hh) = 2/ |/^ri/,ir(|A*(£,.)|V + |A"(£L)r^?2) 



+ Re{{ftnfk)'^\hW*{h)}g3 



(71) 



(72) 



(73) 



where A{Z) and A*''"(£i) denote the propagators of the virtual particles in the direct 
channel and in both crossed channels (labelled t,u, cf. Figl2]). 

The quantities D{ZIr), D^IrCr) can be derived from eqs. fl72|) . fl73l) by the substitu- 
tions 



U ^ Ri. A*'"(£z.) ^ A*'"(£«) 
The kinematical factors are 



'^kj 



^kj ) 



fej,k 



rR 
Jij,k- 



92 
93 



mjmk{pi-Pi+). 



(74) 

(75) 
(76) 
(77) 



We can split the terms depending on the polarization of the neutralino into T-even and 
T-odd contributions: 



s?,(x") = s^/(x;) + sr(x"). 



a,O/~0\ 



(7J 



The T-even contributions depending on the polarisation of the decaying neutralino x' 



Es^(x;) = s^"(zz) + ss"(z4) + ^riz^R) + ^nw + ^TiiRin). (79) 



a,E/ 



are 



9' 



\A{Z)\\R', - L', 



X 



[{ReOlff - (/mO;f 



"L\2-\ „a 



kj 



\o:fm-9^2) 



4g' 



cos^ 6 



w 



2g' 



kj ill "T '^kj 

+ /g7.iA-(4)(20:f + O'l] 

i/^ri/".r[iA"(4)p^?2"-iA*(4)iv] 

5 



+ ^e{(/^*)^(/4)W(4)A«*(4)^3n 



^0) 



^1) 



^2) 
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where the contributions T,'^^ {Ziji), Yfjf^ [inljij are derived from eqs. (IHTI) . (182|) by applying 
the substitutions in eq. flTll) and in addition gf2 3~^ ~9i2 3- 
The kinematical factors are 



9l 

92 



mj{p^ope-){pe+s''), (83) 
mj{p^opf+){pe-s''), (84) 
£/3 = mk[{p^ope+){pi-s'') - {p<iopi-){pi+s'')]. (85) 

The T-odd contributions depending on the polarisation of the decaying neutralino Xj 



E^°(x;) = SS^(ZZ) + S^"(Z£i) + S^'^(Z£^) + S^'^(£i4) + S^*^ (£«£«). (86) 



a,0. 



are 



TfiZZ) 



cos"' ©VK 



|A(Z)p(L^ 2i?e(0,,^)Jm(0,,^)z^7: 

'L cL* L* \t* ( n \ 



cos^ G 



L^i?e<^ A(Z) 



9l 



(88) 

2g'Re[{f^r{fi\?^\h)A-*{h)9t}, (89) 

where the contributions Yfj^ [Zln), Yfj^ {lulji) are derived from eqs. flHTl) . (182|) by applying 
the substitutions in eq. flTll) . The kinematical factor is 

(90) 



^4 



i'm'k^lJ.upcTS"'^p'^op'^-Pl+ 



B.3 Top decay t W+b 

We provide analytical expressions for the 2-body decay of the top quark into a ly-boson 
and the final-state bottom quark: 



2 4 



m 



w 



The spin- dependent contribution is T-even and reads: 



9 r / b/ \ "^t ~ ^ 



-mt{{s\t)pt,) + 



m 



''^'s\t)pw)}. 



(91) 



(92) 



w 



C Kinematics 



C.l Phase Space 

The complete cross section for the process can be decomposed into the production cross 
section and the branching ratios of the subsequent decays: 



daTotai = da{gg iiii] 



txl) 



A2 A2 
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dT{t W+b) , 



(93) 



where the factors E/mV come from the use of the narrow- width approximation for the 
propagators of the t, X2 ^"^^ ^- This approximation is vahd for (F/m)^ <^ 1, which is 
satisfied for Fj ~ 1.5 GeV [43] and F^ ~ 4 GeV. It is also trivially satisfied in the case of 
T^o ~ 10"'^ where the width is small because only the three-body decay is kinematically 
possible. 
We have: 

rfF(ti ^ txl) = -^Pixlt) d<l>i, (94) 

- xTn = ^^(^2) d<l>^o, (95) 

dT{t^Wn) = -^D{t)d^t, (96) 
where the phase-space factors in the laboratory system are given by: 

1 Ip^oP 

= — — - — —r-, — , , ^ dn^, (97) 

* i27T)^2\Ei\p%\-E^Mcosei\ *' ^ ^ 



X2 ' X2 ' 



d$^o = — 1-- Ei-dEi-dQ,+ dQ,-, (98) 

^2 8(27r)5 ||pj^o|cos^^;+ — E^o — Ei+ — Ei-cosa\ 

d<!>t = tA^T^ 7^ ^dn,. (99) 

{27ry 2||pt|cos6';, - Ew - Eb\ 

There is a subtlety in the phase-space calculation, namely that there can be two solutions 
for p^o. If |p^o| < po where po = A2 (m?, m|o, mj)/2m^o, then the decay angle, = 

Z(p£, pt), is unconstrained and there is only one solution. However, if p^o > po, then the 
angle is constrained by sin6'™"^ = Po/|Px§| and there are two physical solutions 

(m? + m|o - ml)\pi\ cos 6^ ± Ei^X{m\, m'io^mj) - 4|pif m|o (1 - cos^ Oi) 
" 2|p,f (l-cos2 ei) + 2m} '^^^^^ 

For the region of phase space where two solutions exist the cross section becomes a 
summation of the solutions for each of the subsequent decay chains. 



C.2 Integration limits 

When evaluating the phase-space integral at the parton level, kinematical limits need to 
be determined on some of the variables and these are listed below. 

If |p^o| < Pq, where Pq = A2(m?,m|o,m^)/2m^o, there are two solutions for p^o, 

eg. (11 00 p . and the decay angle of the t is constrained by 
A5(m?,m|o,m2) 

sin^,- < . (101) 

2|pt|m^o 
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The three-body-decay phase space of the X2 ^^^so has hmits: 

Ee- < (102) 

cosOe- < ^. (103) 
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